Aligned crystallite powder of NdFeAsO0.86F0.14: magnetic hysteresis and penetration depth 
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We report the basal-plane critical current and superfluid density of magnetically aligned NdFeAsOo.geFo.M 
powder. This sample has individual crystallite grains permanently oriented with their c axis along the external 
measuring field. Magnetic irreversibilities at high field suggest strong flux pinning of basal-plane critical cur- 
rents, with monotonic field dependence and no evidence of the "fishtail" effect. The small particles provide a 
sensitive indicator of dc flux penetration, and allow analysis of the temperature dependence of afe— plane London 
penetration depth Xab,i^, which is quadratic at low T . This feature may not necessarily be due to the nodes in the 
gap, but may be rather a sign of a strong pair-breaking. A quantitative determination of the absolute magnitude 
of Aaft L is hindered by the need for accurate knowledge of the particle size distribution. 



PACS numbers: 

I. INTRODUCTION 

Discovery of a new superconducting family of FeAs-based 
superconductors has generated an explosion of research ac- 
tivitjC. The superconductors in this family were produced by 
doping a parent compound, REFeAsO ("111 1" phase, RE-rare 
earth) with electrons by substitution of fluorine for oxygen, 
or by some oxygen deficiency. The highest Tc obtained this 
way is 55 K in^ SmFeAsOo.gFo.i. Hole doping by substitut- 
ing strontium for rare earth also leads to superconductivity^ 
with a Tc = 25 K. Later, superconductivity in a related mate- 
rial BaFe2As2 (these are termed "122" materials) was discov- 
erecP, where again, superconductivity occurs by doping either 
electrons (e.g. Co on Fe site, Tq k, 22 K) or holes (e.g. K on 
Ba site, Tc ~ 38 K). The parent (undoped) compounds of both 
1111 and 122 materials are antiferromagnetic metals with a 
spin-density wave setting in at around 130 k!^. The Fermi 
surface consists of multiple sheets, electron- and hole-like in 
character^. Of the superconducting properties, most basic are 
the two Ginzburg-Landau length scales: the a/j-plane values 
for the coherence length ^^b is about 2-3 nm (e.g. Ref.^ and 
the magnetic penetration depth Xab ~ 200 nm (e.g. Ref.l^i) 

In this work we concentrate on NdFeAsOo.geFo.M super- 
conductor with a 7c = 48 K, as a member of the iron pnic- 
tides 1111-family. Like other pnictide superconducting com- 
pounds, it too has a layered, anisotropic crystal structure. Re- 
cently, the synthesis of small single crystals of such high 7c, 
1111 iron pnictide superconductors has enabled characteri- 
zation of some anisotropic superconducting propertied"''"! 
It is important to further elucidate properties determined by 
strong basal-plane supercurrents. For such purposes, an alter- 
native to small single crystals (typical mass < 1 /xg""' '^) is 
a large collection of small crystallites (here, total mass '--^ 10 
mg), whose orientation has been made to have a nearly com- 
mon c— axis by magnetic alignment. Such samples actually 
have some advantages over small single crystals, since a sub- 
micron grain size allows higher sensitivity to small changes in 
magnetic penetration depth, while the large number of parti- 
cles provides a strong signal. Therefore, it is possible to con- 
duct meaningful measurements using low-frequency and dc 



techniques. In this article we report on making such a sample 
and our findings from both the high field persistent (critical 
state) currents and the low field screening. 

II. SAMPLE PREPARATION AND CHARACTERIZATION 

Stoichiometric amounts of FeAs, NdFs, Nd203, and Nd 
metal were ground together in a helium glove box, pressed 
into a pellet and sealed in a silica ampoule with 0.25 atmo- 
sphere of Ar gas. The pellet with nominal composition of 
NdFeAsOo.86Fo.i4 was heated to 1200° C for 15 h. Powder 
x-ray diffraction gave the correct phase with about 2-5% of 
NdAs, FeAs, and NdOF impurities. The Tc was 48 K, as de- 
termined by the low-field ac and dc susceptibility. 

A polycrystalline piece of NdFeAsOo.geFo.M was mechan- 
ically ground to a fine powder, such that each particle con- 
stitutes a single grain. We did not use any sieve to remove 
large particles. A ten mg quantity of the as-ground powder 
was mixed with clear epoxy (Hardman 04004, setting time 1 
hour) in a cylindrical gelatin capsule and placed in a magnetic 
field of 6 T, which aligns the powder particles with ab planes 
parallel to the field. A similar phenomenon was noted by H.- 
J. Grafe et al^ and B. C. Chang et al.^ although in those 
cases the rare earth ions were respectively spinless La^+, and 
Sm^+ with a magnetic moment of ^ 0.7 jJ-s, while in our case 
it is Nd^+ (^79 /2 ground state, which has a magnetic moment 
of 3.3 jJ-B as a free ion). Therefore such an alignment cannot 
be ascribed to rare-earth magnetism, but rather results from 
layered, anisotropic crystal structure. We have performed the 
alignment procedure at room temperature in two steps to ob- 
tain c-axis orientation, as illustrated in Fig. [T] In the first step, 
field was applied perpendicular to the capsule axis, which pro- 
duced a sample with c axis randomly distributed in the plane 
perpendicular to the applied field (Fig.[TJa)). After about half 
hour, the capsule was rotated 90° about its long axis (fig.[TJb)) 
and the epoxy was allowed to set. As a result, particles reori- 
ented to align with a common c-axis (Fig. [Tfc)), with a pos- 
sible exception of a very small number of particles, whose c- 
axis was perpendicular to the capsule axis after the first align- 
ment step. The total powder volume was ~ 1.3 mm^ and the 
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FIG. 1 : (color online) Two-step alignment of the powder particles, shown asa — b platelets for illustration only. Magnetic field of 6 T points into 
the page. The first step produces particles with ab plane parallel to external field and their c-axis randomly oriented in a plane, perpendicular 
to the field (a). After 90° rotation (b), the second step aligns c with the sample's long axis (c); (d) 6 — 26 scan shows mostly (OOL) peaks and 
predominant c-axis orientation. Background is due to glass-like epoxy. (e) Rocking curve of (004) peak shows 18° FWHM spread in individual 
particles orientation, (f) Absence of 4-fold structure in ij) scan of (114) reflection indicates no azimuthal ordering of powder particles 



approximate capsule volume of 100 mm^ resulted in a dilute 
< 1.5 vol. % particle composite dispersed in the epoxy ma- 
trix. 

Grain alignment was confirmed using a 4-circle x-ray 




FIG. 2: SEM micrograph of a single particle on a holey carbon TEM 
grid, demonstrating diameter ~ 200nm and approximately equiaxial 
shape 



diffractometer with a Cu Ka source (50 kV, 100 mA), a sagi- 
tally focusing graphite monochromator, soller slits, and a Nal 
scintillation detector AO — 29 scan (Fig.[TJd)) shows that the 
(003) and (004) reflections are enhanced relative to a random 
powder and the rocking curve for the (004) reflection is about 
18° FWHM (Fig.[Tfe)). This orientation spread is acceptable 
for studies of in-plane magnetic properties with field along the 
c axis, but is too large if one wishes to study the behavior in 
parallel field. Only the H \\ c situation is examined in this pa- 
per An azimuthal (0) scan at a 37° tilt (Fig. [TJf)) shows that 
the (114) intensity has a weak 2-fold variation which we at- 
tribute to the sample shape, i.e. there is no indication of an 
azimuthal grain alignment, which would produce at least a 4- 
fold symmetry. 

For a quantitative analysis of magnetic data on such sam- 
ples it is necessary to know the particle size distribution. An 
SEM image of a particle on a holey carbon TEM sample 
grid (Fig. |2|i shows a diameter of approximately 200 nm and 
close to an equiaxed particle shape. More detailed informa- 
tion about particle sizes was obtained using a NICOMP ZLS 
paiticle sizing system'^, which deduces particle size distribu- 
tion from temporal coiTelations in intensity of light scattered 
off the particles diffusing in a liquid of a known viscosity and 
refraction index (we used diffusion pump fluid 704). The 
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Averaging with respect to particle size modifies the rhs of the 
above equation: 
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FIG. 3: Magnetic moment of aligned particles in epoxy for field ap- 
plied II ah plane (filled symbols) and [| c axis (open symbols) at room 
temperature. Particles align with a larger magnetic moment along 
the applied field. Inset: grain size distribution function, obtained by 
light scattering from a sample of the powder 

output, a smooth log-normal distribution S{r) (Fig. [S] inset), 
shows that most of the particles are in sub-micron range. In 
a small applied field, superconducting particles smaller than 
a magnetic penetration depth produce little diamagnetism. 
Rather, the particles in the tail of the distribution, in the range 
of about a micron and above, most affect the measurement. 
Therefore the knowledge of the distribution tail is vital for 
quantitative analysis. 

Following the sample preparation we have confirmed that at 
room temperature the magnetization was higher for field along 
the ab plane, than along the c axis, as shown in fig. |3] The 
ratio of Mab /Mr = 1 .26 at // = 1 T is similar to that obtained 
by Chang et al.. Particles align with a larger moment along 
the field. 



III. MEASUREMENTS 
A. Magnetic hysteresis 

The magnetic hysteresis has been measured in dc and os- 
cillating sample (RSO) modes of the Quantum Design MPMS 
magnetometer. The agreement between the two methods is 
good. The Fig. |4|a) shows the RSO data at several tempera- 
tures. 

The ab-plane critical current density was extracted from the 
magnetic hysteresis measurements below 6 Tesla. Knowledge 
of particle sizes is crucial to quantitative analysis. Spherical 
particles were assumed with size distribution of the Fig. [3] in- 
set. Clem and Kogan^ found that the magnetic moment m 
of a sphere of the radius R and critical current density Jc in a 
Bean critical state is given by 



V 320^ J^S{r)rUr 
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where S{r) is the particle size distribution function of Fig. [3] 
inset. 

The extracted Jc are shown in Fig.|4|b). Apparently, values 
for the in-plane critical current are relatively high. For exam- 
ple, at r = 5 K and applied field Bapp = 6 T the critical current 
density is about 0.3 MA/cm^. It should be mentioned, that the 
magnitude of the Jc is very sensitive to small uncertainties in 
the size distribution. The relative values are, of course, quite 
well defined. 

The magnetic field dependence in Fig. |4fb) exhibits a dis- 
tinct region of power-law behavior, Jc °^ with a nearly 
temperature-independent value for a of about 0.38 except 
near Tc, at 40 K. A weak dependence of a on T has been 
found in high-7c cuprate coatings as welP^. Unlike other re- 
ports in the liter ature on the irreversible properties of iron- 
arsenide j'^ l ^^ l ^^ l we do not observe any "peak" or "fishtail" 
effects in Jc{H)- 

Figure|4|c) shows the temperature dependence of the back- 
ground paramagnetism, which is predominantly due to Nd 
ions. This magnetic background has been estimated from 
the hysteresis loops of Fig. |4|a) as follows: first, a mean 
of the upper and lower hysteresis branch is calculated for 
each loop, as a proxy to the equilibrium magnetization M^.^^ w 
1/2(M++M_). These are shown in FigQa) as solid blue 
lines. Then, a tangent to such center hne, passing through 
origin, is constructed (for T — 5 K this is shown in FigQa) 
by the red dashed line). The slope of this tangent provides 
a good approximation for the background susceptibility Jpara- 
The Curie-Weiss fit to these values, ;i^para + (dashed 
line in Fig. |4|c)) produces a positive 9 = 4.7 K, pointing to 
antiferromagnetic coupling between Nd ions and an effective 
Nd magnetic moment of about 2/Zb, i.e. about 60% of the 
free ion value. For T > Tc the magnetic susceptibility 

decreased with increasing T, but could not be accurately de- 
scribed by a Curie-Weiss dependence. We attribute this, as 
well as the reduction in the inferred Nd^+ magnetic moment 
value compared to that of a free ion to crystal-field effects in 
the aligned array of particles. Knowledge of the background 
paramagnetism will be required below to obtain a London 
penetration depth from the measured ac magnetization. 



B. Magnetic penetration depth 

We measured the low-field magnetization M{T), using an 
ac magnetometer (Quantum Design MPMS) between Tc and 
2 K in zero applied dc field (Fig.|5]). The drive amplitude and 
frequency were Hac = 1 Oe and f —\ kHz respectively. The 
epoxy had a weak, temperature-independent diamagnetism, 
about 0.2% that of the sample's signal at low temperature. 
This background has been subtracted from the data. Knowing 
the particle size distribution S{r), from the measured M{T) 
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FIG. 4; (color online) (a) Magnetic hysteresis of the aligned powder samples at (top to bottom) T = 2, 5, 10, 20, 30, and 40 K. Only RSO 
data are shown. Solid blue lines at the average of every hysteresis loop are an approximation for equilibrium magnetization M^q. The slope 
of the dashed red line (shown only for T = 5 K) serves as an estimate of paramagnetic susceptibility of Nd ions, (b) Critical current density, 
extracted from magnetic hysteresis data using Bean model for spherical particles (Eq.|2| for (top to bottom) T = 2, 5, 10, 20, 30, and 40 K. 
Power law dependence of on field is indicated, (c) The temperature dependence of the background paramagnetism, given (for T = 5 K) by 
the slope of the red dashed line in the panel (a), fits well with the Curie-Weiss law (black dashed line), producing 9 = 4JK and an effective 
Nd moment of 2^b- 



we can extract the magnetic penetration depth X^h, using the 
equation 

JoPPir/K,h)S{r)dr 



4nM = - : 



Ha, 



£r^S{r)dr 



(3) 



where, assuming spherical particles, we have a fraction of 
the grain's volume from which magnetic flux is excluded^ 
P{rlKb) = [1 - (3A„/Vr)coth(r/A„i) +3A2yr2] and the de- 
magnetizing ratio D ~ 1/3. This is a straightforward tech- 
nique, which generally agrees very well with other methodd^^. 
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FIG. 5: ac magnetic moment of the aligned powder. Drive amplitude 
1 Oe, frequency 1 kHz. A-in-phase component m' , A-out-of-phase 
component m". There is no peak in m" because most particles are 
small compared to magnetic penetration depth. 



However, as mentioned above, the relatively few large parti- 
cles produce most of the diamagnetic signal. Unfortunately, 
the large r tail of the size distribution is not known with 
enough accuracy to obtain reliable determination of the ab- 
solute value of the penetration depth. In contrast, the temper- 
ature dependence of X^b is far less sensitive to uncertainties in 
size distribution, as shown below, so that we can obtain valu- 
able results for the variation of Xat with temperature. 

Circles in Fig.|6|a) show the X^i^{T) extracted from the as- 
measured magnetization. There is a maximum near 10 K, be- 
low which the diamagnetism is reduced by the underlying Nd 
ion paramagnetism, as is evident also in Fig.|5] Clearly, the Nd 
magnetic signal must be accounted for. Knowing the param- 
agnetic background from the previous section (Fig.|4]^c)), we 
can do such a correction. We demonstrate here that, while 
it is possible to under-correct for the magnetic background 
and make the measurements appear to be consistent with a 
fully gapped state (with two different gaps on different Fermi 
sheets), a proper treatment of the background points toward 
an unconventional pairing or possibly a strong pair-breaking. 

The Nd ion paramagnetism provides two contributions to 
the overall signal, as discussed by J. R. Cooper^". First, the 
paramagnetic signal comes from the grain's field-penetrated 
part, whose volume fraction is given by 1 ~P{r/X). Follow- 
ing to subtract this paramagnetic signal we must re- 
place P{r/X) in Eq. ^ by P^ff "^^Xp^vS - P)- After 
such correction, the data of Fig. |6|a) show no maximum at 
any temperature (green triangles), but rather flat behavior at 
low temperatures, consistent with a fully-gapped state. On 
closer inspection, no simple i-wave BCS calculation fits the 
corrected data satisfactorily. Instead, as shown in Fig.[6|a) by 
the dashed black line, we obtain a good fit by assuming two 
gaps, with T — Q values 2.15 ksTc and l.lAkgTc, the smaller 
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FIG. 6: (color online) (a) The ^^^{T) extracted from the measured 
susceptibility x{T)- without corrections, -corrected for extra 
signal from field-penetrated part, ^-additionally corrected for mag- 
netic permeability. Solid and dashed lines show quadratic and 2-gap 
BCS fits, respectively(b) The T— dependence of the superconducting 
gaps, used in two-gap BCS fit of . (c) Enlarged low-temperature 
region, showing quadratic temperature dependence of the 



gap producing a quarter of the total superfluid fraction. To 
obtain the gap temperature dependence (Fig. |6|b)), the self- 
consistent BCS equation was solved numerically. As men- 
tioned above, the shape of the X^^{T) curve depends only 
weakly on how well the size distribution is known. For ref- 
erence, a change in particle size distribution that would de- 
crease A^^(O) by about a factor of 100, changes the large and 
small gap magnitudes to 2 and 1.23 of ksTc respectively, and 
the relative weight of the smaller gap increases to about 35%. 
Therefore, based on the correction P Pgff alone it could be 
concluded that this superconductor is fully gapped and is de- 
scribed well by the clean-limit BCS theory with two different 
gaps. Indeed, many authors have arrived at such a conclusion, 
as discussed below. 

There is a second correction, however, that must be made 
in order to fully account for the magnetic background, and 
this one comes from the magnetic permeability in the Lon- 
don equation, mV x V x B = — /ione^B, which describes the 
shielding of the magnetic field in the interior of a supercon- 
ductor^*^. For a magnetic mateiial like NdFeAsOF, jiQ must 
be replaced with /z/io = (1 +47rXpara)Mo (factor of 4n is be- 
cause the latter formula is in SI units, while values of fig.|4|c) 
are calculated in Gaussian units, and there is a difference of 
47Z between SI and Gaussian values of magnetic susceptibil- 
ity). Then A^^ = HQue^ /m must be multiplied by 1 +4;r;ifpara 
to obtain the true London penetration depth A„i l, which is 
the measure of the density of the superconducting conden- 
sate. This is shown by blue diamonds in the Fig. |6|a) and 
(c). At low temperatures these are best fit with a parabola 
1 - (r/55K)2 (solid blue line, ¥ig^ c)). Similar analysis and 
result can be found in Ref.^^, from which we have gleaned the 
idea of this second correction. In that work, the role of Pgff 



is negligible due to large size of the sample, most of which is 
shielded from the magnetic field, while the correction for pL is 
large due to a much lower value of the Cuiie- Weiss Q = 0.2 
K, treated as a fit parameter there, resulting in a much stronger 
T-dependence of Jpara near T = compared with the present 
work where Q = 4.7 K has been obtained independently. 



IV. DISCUSSION 

There have been several reports on the temperature de- 
pendence of the superfluid density in various 1111 oxy- 
pnictides. From experimental work one can note measure- 
ments of H. Luetkens et al. (muons, LaFeAsOi_jF;cj^, R. 
Khasanov et al. (muons, SmFeAsOi_j:F;t^ NdFeAsOi_;tFij^ 
on polycrystals, and L. Malone et al. (rf resonator, 
SmFeAsOi_j;Fj:)^-, k. Hashimoto et al. (microwave cav- 
ity, PrFeAsOi_j;F;cj^, and C. Martin et al. (rf resonator, 
NdFeAsOi_rFx, LaFeAsOi-^F^jSl, done on single crystalline 
samples. With the exception of Ref.l^S there is a common 
agreement that the superconducting state is fully gapped 
and can be accounted for by two gaps on different sheets 
of the Fermi surface. This manifests itself in temperature- 
independent A^^ at low temperatures, followed by convex 
temperature dependence. It is possible that this result would 
be different, had the proper corrections been applied. On 
the theoretical front, L. Benfatto et al.^^ have developed a 
model, where the superconducting properties are determined 
by an interband coupling between one electron- and one 
hole-like bands. The model reproduces the main features 
of A^(0)/A^(r), and agrees with the green dashed curve of 
Fig. |6] quantitatively: Ai(0) = 2Tc and A2(0) = Tc- The ex- 
peiiments of L. Malone et al^^ provide similar estimates of 
two gap magnitudes. C. Martin et a/.ISlalso argue in favor of 
two gaps on different Fermi surface sheets, but find a be- 
havior of Xab{T) in single crystals of both NdFeAsOi-vF.T and 
LaFeAsOi_;tFv (where for the latter the paramagnetic back- 
ground is absent and interpretation is more straightforward). 
Their findings are similar to the present result (blue diamonds. 
Fig. |6|a) and (c)) suggesting a possibility of unconventional 
superconductivity. We would caution that in and of itself, the 
quadratic temperature dependence of the superfluid density 
can not be taken as a sure sign of unconventional supercon- 
ductivity. For example, in Nb/Ni bilayers the superfluid den- 
sity displays a quadratic T-dependence due to pair-breaking 
proximity effect between superconducting Nb and feiTomag- 
netic Ni (Ref.^?). With the proposed pairing state of iron- 
pnictideg^the roles of magnetic and non-magnetic scatterers 
are expected to be interchanged: non-magnetic scattering will 
broaden the peak in the superconducting density of states until 
eventually one or both gaps close (perhaps not everywhere on 
the Fermi surface at once), while magnetic scattering will tend 
to make gaps equal in magnitude and reduce gap anisotropy 
by connecting states from different parts of the Fermi surface 
(Anderson theorem). Nagai et al.^^ have analyzed a case of 
anisotropic gap on the electron-like Fermi surface (around the 
M point of the Brillouin zone) and have shown, in particu- 
lar, that the observed power-law temperature behavior of the 
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spin-lattice relaxation rate rj^'(Ref!^, considered by many 
to be at odds with a fully gapped state, finds a plausible ex- 
planation in that approach: while the density of states is fully 
gapped {N = for |£| < Ai), it varies linearly with energy in 
the interval Ai < |£| < A2, a feature reminiscent of the cZ-wave 
situation in cuprates. We note that a qualitatively similar be- 
havior of A^(e), i.e. a finite sub-gap quasiparticle density of 
states could result from non-magnetic scattering. Overall, the 
question of the pairing state in iron oxypnictides is far from 
resolution at this stage. 

There have been several reports on the critical current den- 
sities of 1 1 1 1 materials, all extracted from magnetic hyst eresis 
of polycrystalline samples by means of the Bean modeP^EH 
The shortcoming of this approach is that the spatial scale of 
supercurrent circulation is often unknown: it can be given 
by an individual grain size (probably the case near 7c and/or 
in high fields), or by the sample size, which is likely at low 
temperatures and fields in well connected materials'^'*. The 
resulting inferred critical current density can differ by two 
orders of magnitude. In a single crystal of the 122 ma- 
terial Bai_TKxFe2As2, the in-plane critical current density 
reaches^^ close to 5 • 10^ A/cm^ at 2 K, T, a value com- 
parable to our result. 

V. CONCLUSIONS 

We have obtained the temperature dependence of the 
plane superfluid density ^^y]^{T) and critical current 



Jc{T,H) as a function of temperature and magnetic field 
for an ensemble of aligned, sub-micron size grains of 
NdFeAsO0.g6F0.14. The temperature dependence of the super- 
fluid density, after proper corrections for background param- 
agnetism, is quadratic, which alone does not enable us to come 
to a conclusion regarding the nature of the pairing state. The 
critical current density is high, which makes this superconduc- 
tor potentially interesting for applications, such as supercon- 
ducting wires. We have not observed any "fishtail" anomalies 
in Jc{H). 
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